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Abstract. Saturation is expected to occur when a high density of partons (mainly gluons)- or 
equivalently strong fields in Quantum Chromodynamics (QCD) - is realized in the weak coupling 
regime. A way to reach saturation is through the high-energy evolution of an extended target probed 
at a fixed hard scale. In this case, the transition to saturation is expected to occur from nonlinear 
perturbative QCD dynamics. We discuss this approach to saturation, which is mathematically 
characterized by the appearance of traveling wave patterns in a suitable kinematical representation. 
A short review on traveling waves in high energy QCD and a first evidence of this phenomenon in 
deep-inelastic proton scattering are presented. 



INTRODUCTION 

One of the most intringuing empirical observations about deep-inelastic scattering 
(DIS) on a proton is the geometric scaling property yj], which states that the y* — p 
total cross-section data can be approximately ploted on a one-dimensional curve 
0T*~/>(<2 2 /Q 2 (Y)) where (Q 2 ,Y) span the virtuality-rapidity kinematical domain of 
DIS. One finds Q 2 (Y) » e XY . 

The question we want to address is whether the observed geometric scaling can be 
explained in terms of QCD evolution equations describing saturation. 

Saturation of parton (mainly gluon) densities at high rapidity may be at the origin of 
the geometric scaling property. Saturation is expected to occur when a high density of 
partons is created in the limited geometry of the target. In the expected scenario, the 
parton wave functions overlap and eventually create a new state of matter, the Color 
Glass Condensate (CGC) [2]. In the case of DIS on a proton, increasing densities are 
obtained by high-energy evolution at fixed virtuality Q 2 , when the corresponding linear 
QCD evolution, named BFKL after Balitsky, Kuraev, Lipatov, and Fadin [3], gives 
rise to an exponential growth of the number of gluons. At high enough energy the 
BFKL formulation gets modified. In the approximation of uncorrelated hard probes, 
one deals with the Balitsky-Kovchegov (BK) equation o], where the corrections due 
to the high density of partons are expressed in terms of a non-linear damping term. We 
shall investigate whether and how geometric scaling could be due to the nonlinear QCD 
evolution of the gluon distribution in the target. 




FIGURE 1. Schematic picture of the transition to saturation with rapidity. In a region Y ~ Y\ : the 
projectile of size ~ l/Q is able to probes the number of partons in exponential growth; when Y ~ Y2 : 
the probe counts partons by groups, leading to a nonlinear damping evolution factor in the probed parton 
density. For Y > Y2 : A new phase of partonic matter (CGC) may appear through parton correlations. 

GEOMETRIC SCALING AS TRAVELING WAVES 

Let us start from the BK equation expressed in momentum space (using the impact- 
parameter independent formulation) 

d Y ^V = ax (-d L ) JV - a ^ 2 

where jV(J,L = 21ogfc) is related to the tranverse-momentum k distribution of gluons 
in the target. The (leading order) BFKL kernel is 

x (-d L )=2 ¥ (l)- ¥ (-d L )- ¥ (l+d L ) . 

The coupling constant will be considered either constant, or running like a = 1 / bL. 

From a mathematical point of view, the BK equation can be related [(J to an "univer- 
sality class" of nonlinear equations for which characteristic properties can be rigorously 
derived. "Universality" means here that the solutions have general characteristics which 
are essentially independent from initial conditions. In the case of the BK equation (and 
also for some extensions beyond BK) this generic feature is the formation of traveling 
wave patterns when the energy increases. As we shall see, traveling waves are intimately 
related to geometric scaling. Let us define these traveling waves and give a qualitative 
explanation of their formation during the rapidity evolution. 

Consider first the "diffusive approximation" of the BFKL kernel 
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FIGURE 2. Typical traveling wave pattern. The wave front (in QCD: amplitude) connecting the regions 
u = 1 (in QCD: strong absorption) and u = (in QCD: transparency) travels from the left to the right as t 
(in QCD: Y) increases. 



By a suitable redefinition of the function and variables one can map the BK equation in 
the diffusive approximation onto the Fisher Kolmogorov-Petrovski-Piskounov (F-KPP) 
equation ED describing a reaction-diffusion process in space and time: 

d t u(t,x) — d 2 u(t,x) + u(t,x) —u 2 (t,x) , 

In reaction-diffusion language, d 2 u(t,x) is the diffusion term, u(t,x) is responsible for 
the exponential growth and u 2 {t,x), the damping term. 

The "dictionnary" between F-KPP and BK can be written as follows 

Time t — > Y 

Space x -> L+\a%" (\)Y 

Wave Front u(x — ct) — > ^V(L — vY) 

Traveling Waves — > Geometric Scaling. 

The key feature of the solutions of the F-KPP equation (and of equations belonging 
to the same universality class) is the formation of traveling wave patterns, see Fig.©. It 
can be rigorously proven that, at large times (=>■ large rapidities), the solutions verify 

u(t,x) > u(x-m(t)) jY{Y,L) > jy(L-L s (Y)) 

t — > oo Y — > °° 

with 



L S (Y) = logQ^(F) = a X '(f)Y - Alog7 - ±, J + . 



7 is the implicit solution of the characteristic equation x(f)/Y = X'(7) which leads to 
the asymptotic wave speed v = ocx'if)- It plays the role of a critical QCD anomalous 



FIGURE 3. "Pulled vs. Pushed fronts". The function u(t,x) is represented for the three different classes 
of initial conditions. Top: "Pushed fronts," the wave front keeps the initial front profile with a non- 
universal speed v_; Middle: "Pulled = Pushed"; Bottom: "Pulled front:" both the wave front and the 
speed acquire universal values at large t (~ Y) in an intermediate kinematical region called the "wave 
interior". 



dimension. One immediately recognizes the equivalence between geometric scaling and 
traveling waves solutions of BK. Moreover, quite a detailed information can be obtained 
from the sole knowledge of the linear kernel, as for the three first terms [6] of the 
asymptotic expansion of the saturation scale Q% (Y). 

ASYMPTOTIC TRAVELING WAVES 

In fact the results using the diffusive approximation of BK and its mapping to the F-KPP 
equation are more general. Instead of entering here in refined (and technically useful) 
mathematical arguments [8], let us qualitatively illustrate the reason why traveling wave 
patterns generically develop during the rapidity evolution from a given initial condition, 
using the concepts of "pulled fronts". 

Universality features are obtained in the "pulled" front situation, see the bottom plot 
of Fig.®. In this case, one starts with initial conditions corresponding to an anomalous 
dimension 70 > f, where 7 = .6275... is for the (leading order) BFKL kernel [6]. The 
formation of the traveling wave comes from the competition between the exponential 
growth at small values of the solution u, where the linear term dominates the evolution, 
with the nonlinear damping exercised at larger values of u. Constrained by these two 
opposite trends, the solution is "forced" to adopt an universal behaviour. Fortunately 
enough, the perturbative QCD initial conditions, through the transparency property 
(70 = 1 > 7), fall into the universal "pulled" front regime |0]. 
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FIGURE 4. The "Reduced" Front profile. 



In Fig.© the "reduced front" (k 2 / 'Q s (Y) 2 y jY \L,Y) , obtained by factorizing out the 
wave propagation, is represented for fixed (left) and running (right) coupling constant. 
It shows the scaling straight line, which envelopes the curves obtained for increasing 
rapidity. Hence geometric scaling violations are also predicted due to diffusion. One 
may notice that the diffusive approach to scaling is slower for running coupling, corre- 
sponding to "anomalous diffusion" in l/? 1 / 3 ~ 1 /y 1//6 (for running coupling t ~ F 1 / 2 ). 
Normal diffusion l/? 1 / 2 ~ 1/y 1 / 2 is characteristic of the F-KPP universality class or 
BK with fixed coupling. 



PARAMETRIC TRAVELING WAVES 

The abovementionned mathematical properties have the major interest of revealing the 
link between geomeric scaling and the nonlinear dynamics of QCD evolution equations 
at high energy. However the path from phenomenology to theory is not yet accom- 
plished. On the phenomenological ground, the simulations of BK equation solutions 
12] have shown that pre- asymptotic s effects may be important and endenger geometri- 
cal scaling predictions from QCD in the physical region where it is observed. On the 
theory side, it has been recently realized [10] that the effect of parton fluctuations and 
correlations, especially in the dilute (transparency) region which is the driving force of 
the pulled front, may be responsible for a breaking of geometric scaling. In fact, trav- 
eling wave patterns still exist lITlll but their event-by-event fluctuations tends to break 
geometric scaling in the average. 

Despite these problems, the numerical simulations [9] also show that patterns behav- 
ing like traveling waves exist in a preasymptotic region in energy. Moreover the mean 
field BK equation seems to be still valid in this region. In order to explore this possibil- 
ity and understand its origin, we proposed id to attack the problem of traveling wave 
solutions of the BK equation in a new way. 

The initial mathematical idea IIOl is to assume and directly insert a traveling wave 
solution u(x,t) — > U(s = x/c — t) into the F-KPP equation, or equivalently to impose 
a geometric scaling form to a solution of the BK equation. Keeping for simplicity the 
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FIGURE 5. The amplitude JV{L,Y) as a fonction of L for different values of F = 3,4, .., 10. Full lines: 
MRST gluon parametrization. Dashed lines: Parametric QCD traveling waves. 



space-time language and the FKPP equation (the dictionnary to BK being still easy), one 
gets 

| dU | 1 d 2 U _ Q 

ds c 2 (ds) 2 

taking into account the previous study showing CD that the wave speed c > 2, we obtain 
an iterative solution 

U(s) = U Q +iu 1 +l t -l-U 2p 



p>2 



obeying an exactly solvable 11211 hierarchy of equations for the U2 p , s. The first one only 
(for Uo) is nonlinear and can be exactly solved, the other boiling down to a rather simple 
linear algebra. This mathematical result translates into a parametric form of a geometric 
scale invariant gluon amplitude jV{L — L S (Y)) which can be used in phenomenology. 
j¥ depends only on few parameters determined by the linear kernel. One obtains 
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This solution (valid also for the running case with an appropriate mapping 11211 ') is 
mathematically valid [8] in the "wave interior", see Fig.©. Remarkably, we observed 
that the corresponding region has a large overlap with the physical region where geo- 
metrical scaling is valid at HERA. In order to check this opportunity, we determined the 
phenomenological observable related to jV{L, Y) from a MRST parametrization of F2 
data. The obtained kernel parameters are not far but different from those expected from 
the leading order BFKL kernel one and thus could require higher order contributions 
JI2I1 . Interestingly, higher order contributions plays an important role in a different but 
related problem: the existence of "backward" traveling waves appearing [14] from the 
energy evolution starting from an inital high density (and not initial low density as in 
DIS) state, which could be the case in heavy-ion reactions. 



CONCLUSION 



Traveling wave patterns have been shown to result from nonlinear QCD equations. The 
phenomenologically observed geometric scaling in DIS on a proton at large energy is 
consistent with the traveling wave patterns. It remains to be found whether one can relate 
the geometric scaling curve with the QCD kernel and more generally to a complete 
solution of high-energy QCD, including correlation and fluctuation contributions. 
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